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Abstract 

A lattice formulation of the 0(1, 2)/0(2) x sigma model is developed, based on 
the continuum theory presented in the preceding paper. Special attention is given to 
choosing a lattice action (the "geodesic" action) that is appropriate for fields having 
noncompact curved configuration spaces. A consistent continuum limit of the model 
exists only if the renormalized scale constant fi^. vanishes for some value of the bare 
scale constant /3. The geodesic action has a special form that allows direct access to 
the small-/? limit. In this limit half of the degrees of freedom can be integrated out 
exactly. The remaining degrees of freedom are those of a compact model having a 
/^-independent action which is noteworthy in being unbounded from below yet yielding 
integrable averages. Both the exact action and the /3-independent action are used to 
obtain /3/j from Monte Carlo computations of field-field averages (2-point functions) 
and current-current averages. Many consistency cross-checks are performed. It is found 
that there is no value of /3 for which vanishes. This means that as the lattice cutoff 
is removed the theory becomes that of a pair of massless free fields. Because these 
fields have neither the geometry nor the symmetries of the original model we conclude 
that the 0(l,2)/0(2) x model has no continuum limit. 
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1 Introduction 



In the preceding paper (denoted here by I) we have outhned the main theoretical features 
that would characterize the quantized 0(1, 2)/0(2) x Z2 sigma model if it had a continuum 
limit in four dimensions. In this paper we present evidence that the model does not, in fact, 
have such a limit. The equations of paper I are therefore purely formal and serve only to 
motivate our procedures in the lattice simulation. 

Because lattice simulations of noncompact sigma models have not been performed previ- 
ously, and because there is evidence that the concept of "universality" of the scaling critical 
point may not apply to noncompact models, we take special care to develop a lattice action 
that is appropriate for the 0(l,2)/0(2) x Z2 model. The chief goal of our investigations, 
which is to determine whether or not there exists an interacting continuum limit of the 
model, can be reached by computing the dimcnsionless renormalized coupling constant /Sij 
(eq. (5.12)). An acceptable continuum limit exits only when /9r is zero. From our numer- 
ical simulations we determine (i^ both from two-point functions, via the effective action 
formalism, and from current algebra. 

In addition to performing direct numerical simulations of the 0(l,2)/0(2) x model 
we also develop a useful small-/9 or "strong coupling" expansion, which permits half of the 
degrees of freedom to be integrated out exactly. The remaining degrees of freedom form a 
relatively simple (although somewhat unusual) compact model which is independent of the 
bare coupling constant. Simulations of this model provide yet another way to study the 
continuum limit of the full sigma model. 

Our numerical results from Monte Carlo simulations of both the full model and its small-/? 
limit indicate that is never zero and hence that there is no interacting continuum limit of 
the theory. Our results are consistent with the existence of a noninteracting continuum limit, 
i.e., a continuum limit that is a free field theory. This result lends support to the popular 
conjecture that all scalar field theories in four dimensions are "trivial" (in the technical 
sense of triviality). However, since this continuum limit has neither the symmetries nor the 
interactions of the original classical model we conclude that a consistent quantum theory of 
the 0(1, 2)/0(2) X Z2 model does not exist. 
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2 The massless scalar field. Lattice action and lattice 
Fourier transforms. 



To get our bearings on how to introduce the lattice we look first at a completely soluble 
system: the free massless scalar field. This system is described by a single real scalar field 
0, and its Euclidean action is 

S=W U,^(t),,d'x. (2.1) 



(We here drop the subscript E introduced in I to distinguish the Euclidean action from the 
Minkowski one.) The field is chosen to be dimensionless and the scale factor is made 
explicit in order to maintain as close a resemblance as possible to the sigma models. The 
lattice is most simply chosen to be a cubical one having N'^ sites and periodic boundary 
conditions (4-torus). On such a lattice the action takes the differenced form 

S = \^^a? ^ [0(one end of link) — 0(other end of link)]^ (2.2) 

links 

where a is the lattice spacing and the links are between adjacent sites in the four principal 
lattice directions, with periodicity conditions taken into account. 
The size, or periodicity distance, of the lattice is 

L = Na, (2.3) 

but to the computer this length is irrelevant. The only adjustable constants in the lattice 
version of the theory are the dimensionless quantities and (3, where f5 is the coefficient in 

13 = fi'^a^ . (2.4) 

If /i is to be a finite mass scale then (3 must vanish as the lattice spacing goes to zero. 

The system ( p.2| ) is best studied by introducing the lattice Fourier transform. Denote 
by (paiBjS the value of the field at the lattice site having coordinates aa, (3a, 7a, Sa, with 
the origin of coordinates taken at a corner of the periodic cube, the coordinate axes being 
oriented along the principal lattice directions, and a, (3, 7, 6 being integers ranging from 
to — 1. The lattice Fourier transform is defined by 

"27ri 



a,/3,7,(5 



j^-{ka + 1(3 + m7 + n5) 



(2.5) 



where fc, /, m, n are integers ranging from —N/2 + 1 to N/2 {N being here assumed to be 
even). The original periodic field may be regained by taking the inverse Fourier transform 
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k,l,m,n 



27ri 
'1^ 



{ka + 113 + mj + n6) 



(2.6) 



It is well known that imposition of periodicity in Euclidean time corresponds to studying 
the field (f) in a thermal state at temperature T = 1/kL, where k is Boltzmann's constant. 
However, since the cutoff energy arising from the lattice discretization of the system cor- 
responds to a temperature of order 1/ka, which is times as big, T may effectively be 
regarded as vanishing, the thermal state being indistinguishable from the ground state when 
is large. 

In order to express the lattice action ( |2.2| ) in terms of the Fourier transform (pkimn, it 
is convenient to introduce some abbreviations. The ordered set of integers {a, (3, 7, 6) will 
be denoted simply by a and the ordered set {k,l,m,n) by k. The following additional 
abbreviations are obvious: 



k ■ a = ka + ip + rri'y + n6 , 
a + a' = {a + a',p + p',-f + y,S + 5') 
— k = {—kj—lj—rrij—n), etc. 



(2.7) 
(2.8) 
(2.9) 



Now let the symbol A (for links) range over the eight values (±1,0,0,0), (0, ±1,0,0), 
(0,0, ±1,0), (0,0,0, ±1) and regard all integers as defined modulo N. Then, taking note 
of the double counting of links, rewrite expression ( |2.2| ) in the form 

CK.A 



a,A,k,k' 



X 



^-{2Tvi/N)k-a. _ ^-{2iTi/N)k-{a+X) 



g-(27ri/Af)k'-Q: _ g-(27ri/Af)k'-(Q;+A) 



where 



A,k 



2 — g{2vri/^)k-A _ g-(27rj/Af)k-A 



Kik) 
K'{k) 



( irk nl Tim 7in\ 
= (^2sin—, 2 sin-, 2 sin— ,2 sin— J , 

= K(k)-K{k). 



(2.10) 

(2.11) 
(2.12) 
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The last line of ( |2.1CI| ) follows from the easily verified identity 

K\k) (2.13) 



^ 1 _ g(27rj/Af)k-A ^ 

^ A 



^-(2TTi/N)k-X 



X 

and the fact that (/)_k = 0k*- 

Note that although the physical continuum momentum corresponding to the integers k, 
I, m, n is K{k.)/a where -ft'(k) = (27r/A^)k, it gets replaced on the lattice by K{k)/a, which 
equals K(k)/a only in the limit N oo. The sine functions appearing in are lattice 

artifacts stemming from the replacement of derivatives by differences. 

With the action recast in the form (|2.10|) it is an easy matter to compute 2-point functions, 
with averages defined as in eq. (1.10.5) but adapted for the lattice. For the free field the 
integrals in (1.10.5) are Gaussian. Since the transformation from the variables (pcx to the 
variables 0k effectively diagonalizes the Gaussian matrix, and since this transformation is 
unitary with unit Jacobian, one may write down the 2-point function by inspection: 

5k,-k' 



The general ra-point functions are equally easy to obtain. They vanish when n is odd and 
are expressible as sums of permuted products of 2-point functions when n is even, which is 
another way of saying that all higher correlation functions vanish for Gaussian probability 
distributions. 

Actual machine computations of the 2-point functions are found to agree very precisely 
with eq. ( p.l4| ) for all values of k and k' except (0,0,0,0). The reason (0,0,0,0) is special 
(apart from the fact that expression ( p.l4| ) becomes singular at this value) is that 0o repre- 
sents the constant component of 0, being A^^ times the average of over the lattice. The 
Euclidean action attains its lower bound, namely zero, whenever is a constant. Moreover, 
the addition of a constant to leaves the action unchanged, and hence the integrals over 
the variable 0o diverge. On the computer this divergent degree of integration freedom is 
eliminated by adding an appropriate constant to at the end of each Monte Carlo "sweep" 
of the lattice, in order to keep 0o equal to zero. This has the consequence that the Monte 
Carlo average <^0o") vanishes for all positive integers n. 0o could, of course, be held fixed 
at any other value. Zero is chosen merely for convenience. This fixing is equivalent to 
fixing the vacuum state, the vacuum being here degenerate just as it is in the case of the 
0(l,2)/0(2) X Z2 model 

If one takes the lattice Fourier transform of ( p. 141 ) (remembering 0o = 0) to obtain the 
2-point function directly in ordinary space, one finds that when is large and {ex. — ck')^ is 
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small compared to A^^ the 2-point function accurately approximates the continuum Green's 
function: 

{(l>aM--rYwA < (a - a')' « iV^ (2.15) 

Note in particular that this function is spherically symmetric. For larger values of (a — o:')^, 
to be sure, the cubical toroidal lattice structure begins to show itself, but it is nontrivial 
that ( p.l4|) yields ( |2.15| ) when (a — cx'Y is small despite the fact that -ft'^(k) itself is not 



spherically symmetric. A remnant of 0(4) symmetry remains in the lattice simulation. We 
shall see later that the momentum-space 2-point function of the 0(1, 2)/0(2) x Z2 model can 
be accurately represented by an expression of the form ( |2.14| ) multiplied by a nonsingular 
slowly varying function of K'^{k). We regard this as a lattice confirmation of the 0(4) 
symmetry, or original Lorentz symmetry, of the continuum model. 

3 A difficulty. The 0(1, 1)7^2 x Z2 model. 

The success achieved in simulating the massless free field on the computer prompts one to 
attempt an analogous simulation of the 0(1, 2)/0(2) x Z2 model, based on the similarity of 
expression ( p.l| ) to the part of the action (1.2.5) that does not involve the Lagrange multiplier. 
This similarity suggests that the 0(l,2)/0(2) x Z2 lattice action be chosen in the form 

^ ~ k^'^ Vab 0"(one end of link) — 0" (other end of link) 

links 

(one end of link) — 0^ (other end of link) , (3.1) 

with the (j)"' subjected to the constraint 

r/a6 0V = -l (3.2) 

(cf. (1.2.8)). In the case of the 0(3)/0(2) sigma model, in which 0° is replaced by 0^, r]ab by 
6ab, and the constraint by 0"0" = 1, this is exactly the choice of lattice action that is usually 
made. The reason for this is that if the number of dimensions is reduced from four to three 
the theory becomes formally identical to the statistical theory of the Heisenberg ferromagnet, 
which is interesting to study in its own right: The action becomes proportional to the spin- 
interaction Hamiltonian, /3 becomes the coupling constant divided by temperature, and the 
functional integral of e"'^'''^' becomes the partition function. 

There is, however, a difficulty with the choice ( |3.1| ). To show this we digress for a 
moment by looking at a simpler model, the 0(1, 1)/Z2 x Z2 model, which is obtained from 
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the 0(l,2)/0(2) X Z2 model by omitting the field <^^. Equations (1.2.5) and ( ^72|) continue 
to hold, but r]ah is now given by {j]ah) = diag(— 1, 1). The configuration space is topologically 
IR, instead of IR^, and can be parameterized by a single variable s, where 

0° = cosh s 
(jr" = sinh s. 

In terms of s the Euclidean action becomes 

S = \^i^ j s^^s,f,d^x (3.4) 

which is that of the massless free field! 

The use of two variables, 0° and 0^, to describe the massless free field allows us to view 
it in a new light, and to ask questions that would not have occurred to us if we had confined 
our attention to only one variable, questions that are relevant also for the 0(1, 2)/0(2) x Z2 
model. We note first of all that there is a conserved current, which may be read off directly 
from eq. (1.10.20): 

ioiM = /Lt^ (00 - 01 0O,m) 

= /x^ [— cosh s (sinh s)^^ + sinh s (cosh s)^^] 

= . (3.5) 
The Fourier transform of this current is simply 

Joi/.(p) = ^/^V/.s(p) , (3.6) 
which, in view of the well-known free-particle propagator 

ms(p')) = (3.7) 
leads immediately to the following analog of eq. (1.10.18): 

(Joim(p)Joi.(p')) = f^'^^Sip + P') . (3.8) 

The bare and renormalized Planck masses are seen to be identical for this model, as is to be 
expected. 

As shown in paper I the renormalized Planck mass can also in principle be obtained from 
any 2-point function, for example from {(f)^{x)(j)^{x')). This, however, requires the evaluation 
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of {4>^), and this, in turn, requires introduction of a cutoff such as the lattice provides. Let 
us assume that the correct lattice action for the 0(1, l)/Z2 x Z2 model is expression (|2.2| ) 
with replaced by s. To obtain (0") it is convenient first to compute (s^), using 

,^ = iV-2^'5ke-(2-/^)^«, (3.9) 

k 

where the prime on the summation sign indicates that the term with k = is to be omitted, 
on the assumption that so is held equal to zero. Note that holding so = by adding a 
constant to the free scalar field at the end of each Monte Carlo sweep is equivalent to applying 
an 0(1, 1) "boost" to the field at the end of each sweep. This is completely analogous to the 
procedure that must be adopted in the case of the 0(l,2)/0(2) x Z2 model in order to fix 
the vacuum state. In that case an 0(1, 2) boost is applied at the end of each sweep in order 
to maintain the lattice averages of (p^ and 0^ equal to zero. 

Using eq. (|2.14|) , with replaced by s, and invoking the displacement invariance of the 
theory, one gets 



a k,k',a 



N-'Y.'i~'^^-^) = FiN)//3 (3.10) 



k 



where 



F{N) = AT-^^ ir-2(k). (3.11) 

k 



In the Appendix the following easy estimate of F{oo) is derived: 

F(oo) ^ .15. (3.12) 

This estimate has been confirmed by explicit summation on the computer. F{N) is found to 
converge rapidly to the value .1549. . . for greater than 4. This is an example of a stability 
property that one frequently encounters in lattice simulations: Important quantities often 
(although not always!) become independent of A^ as A^ is increased. 
Because a free field obeys Gaussian statistics one has also 

whence 



^ ) = S(^ / ' (^^"'^)=0' (3-13) 



(cosh = V At (^'") = E A (^')" = e''^'^^^'^ (3-14) 



and 
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= (sinhs) = (3.15) 
(cf. (1.7.1)). 

Suppose now that instead of proceeding as above one chooses for the lattice action ex- 
pression ( ^.11 ) with a, 6 G {0, 1}. In terms of s this action becomes 

S = 2(3'^ sinh^ | [s(one end of hnk) — s(other end of hnk)] . (3.16) 

links 

When the hnk differences are smaU, expressions ( ^.l(j| ) and ( |2.2| ) (with replaced by s) are 
approximately equal. However, in the functional integral the link differences do not stay small 
in the continuum limit a ^ 0, (3 —* 0. Consequently the two actions lead to very different 
results on the computer. This may be seen as follows: Let As be the difference between the 
values of s at opposite ends of one of the links in the sum (|3.16|) . The contribution that this 
link makes to the sum is 2/?sinh^ long as this contribution is of order 1 or less for 

every link, the field will contribute significantly to the functional integrals in (1.10.5) (with 
replaced by s). The dominant fields are therefore characterized by the condition 

As ~ 2sinh-^^= — ^-ln/3. (3.17) 



Since the cross-link differences are of order — In (3 for small (3, one may expect ( |s| ) to be of 
similar order: 

(|s|) ~ -ln/3 (3.18) 

whence 

_ e-Cin/3^^-C (3.19) 



for some (. Actual machine computations, using the action ( p.l|) , show that ( is approxi- 
mately equal to .57. 

The difference between expressions ( |3.14| ) and ( |3.19| ) for (0*^) is striking. According to 



( p.l9| ) (0°) exhibits a weak singularity as /5 — > 0; according to (|3.14] ) it exhibits an essential 
singularity! Another difference shows up in the computation of the "lattice specific heat," 
which is defined by 

CL = (^{e') - (e)'^ , E^S/(3. (3.20) 

Using the action (2.2) (with replaced by s) one obtains trivially 
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cl = const, (3.21) 

On the other hand, computations based on ( |3.1| ) are found to yield 

Cl = const. x/?-i-s^(3) ^ ^3 22) 



Both ( |3.21| ) and ( |3.22| ) show that /3 = is the "critical point" of the theory. But the critical 



exponents differ! Q 

These results have immediate implications for the 0(1, 2)/0(2) XZ2 model. Obviously one 
cannot reach trustworthy conclusions about any theory, on the basis of lattice simulations, 
unless one can trust the lattice action. 

4 The geodesic lattice action 

The sensitivity of the above results to the choice of lattice action stems from the fact that the 
fields in noncompact models undergo huge fluctuations as /5 — > 0. For the 0(1, 2)/0(2) x Z2 
model we should like to use an action that is as close as possible to the action ( p.2|) (with 
replaced by s) for the 0(1, 1)/Z2 x Z2 model. We first note that the link differences in 



eq. ( p.2| ) are taken directly in configuration space (the real line) whereas those in eq. ([O 
are taken in the embedding space. In our view one should stay in configuration space. But 
unless the configuration space is flat, as it is for the 0(1, l)/^2 x ^2 model, this criterion is 
insufficient. One must also deal with the problem of what "preferred" coordinates to choose, 
if any. 

That this problem may be a critical one can be seen by noting that the dominant cross- 
link differences to which the action ( p. 21) gives rise satisfy 



A0~^ = — . (4.1) 

Even if one imagines the continuum fields to be C*°° interpolations of the lattice fields, one 
sees that in the functional integrals these fields become wildly discontinuous in the limit 
a — ^> 0. Their derivatives are even more singular, diverging nearly everywhere like 1/a^ as 
a ^ 0. 

It is easy to see that in d dimensions the derivatives of the fields diverge nearly every- 
where like a"'^/^ as a — > 0. Thus even when d = 1, i.e., in the case of ordinary quantum 
mechanics, the dominant trajectories are nondifferentiable. This is a well known feature of 



^It is also worth noting that ^ cx3 at the critical point already when N is finite. With compact models 
Cl diverges at the critical point only in the "thermodynamic limit" N —f 00. 
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path integration, but it has consequences that are sometimes forgotten. Classically, a system 
may be described by any convenient variables, and in a transformation from one set of vari- 
ables to another, derivatives of one set are obtained from those of the other by application 
of the chain rule. The chain rule is valid, however, only when the derivatives actually exist. 
In the case of path integration in ordinary quantum mechanics it has been known for years 
that the transformation laws for lattice derivatives do not involve merely simple Jacobians. 
Extra terms have to be added to the naive finite differences suggested by the continuum 
classical action. In the case of nonlinear sigma models in four dimensions, in which there 
are no obviously preferred field variables, we do not have any already-worked-out theorems 
to help us in choosing an appropriate lattice action. The quantum theory of these models 
simply does not yet exist as a coherent discipline. We are in the position of having to define 
the theory by choosing a specific lattice action a priori. 
We have chosen to replace the lattice action by 

S=1(3J2A'^ (0(one end of link), 0(other end of link)) (4.2) 

links 

where A is the geodetic distance (1.3.8) in the configuration space {not in the embedding 
space) between the field values at the two ends of the link. This choice has four attractive 
features: 

1. It can be expressed in terms of any convenient field variables, being independent of 
which ones are used. 

2. It reduces to (|2.2| ) (with replaced by s) in the case of the 0{1,1)/Z2 x Z2 model. 

3. It is invariant under the actions of the global group, whether 0(1,1) or 0(1,2) or 
whatever. 

4. When d = 1 (ordinary quantum mechanics) it defines, in the continuum limit, a system 
whose wave function ip satisfies the (almost) obviously correct Schrodinger equation 

td^/dt=^-^n + ^B}i^, (4.3) 

where □ is the covariant Laplacian operator on the configuration space and R is the 
curvature scalar. 



When d = 1 the choice (|4.2|) is also obtained if one simply interpolates classical solutions 



of the dynamical equations between the lattice sites and inserts the interpolated into the 
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continuum action. For c? > 1 a strict analog of this procedure would require interpolation of 
solutions of the classical field equations across lattice cells, but we believe that construction of 
such an analog would be unnecessarily complicated and would yield results not significantly 
different from those obtained with ( [4 .21) . We note also that historical controversies over the 
term in R in (|4.3|) are irrelevant here because, for sigma models, i? is a constant and affects 
only the energy zero point. 

The action ( [4.2| ) can also be used for the 0(3)/0(2) model. But in this case the cross-link 
differences A(f)°- can never be larger than 2, and simulations we have run show that both ( |4.2| ) 
and the conventional action yield results that are qualitatively similar. Both lead to phase 
transitions as P is decreased: At a certain critical value Pc, of order unity, the average (0^), 
which plays the role of magnetization, drops abruptly to small values (which tend to zero 
as ^ oo). The precise value of (3c depends on which lattice action is used, but otherwise 
there is little difference. 

There is a quantity that plays the role of "magnetization" for the 0(l,2)/0(2) x Z2 
model, namely {(j)^) ^. Both ( ^TTD and (|4.2| ) cause (0°) ^ to drop to zero as /3 — > 0, but 
the nature of the dropoff differs profoundly from one action to the other. The reason the 
0(3)/0(2) model has a phase transition is that its configuration space is compact. (0^) drops 
to zero as soon as the functional integral begins to probe the whole configuration space. The 
configuration space of the 0(1, 2)/0(2) x Z2 model, being noncompact, is never fully probed 
as long as (3 is finite. There is always some residual "magnetization" (0°) ^ and there is no 
phase transition. 

These facts have certain implications for the concept of "universality." It is generally 
believed that, insofar as critical exponents near phase transitions are concerned, dimensions 
and symmetries, not details, are what are important. Our insistence on choosing the lattice 
action carefully may seem superfluous. In the case of the 0(3)/0(2) model the choice of 
lattice action probably is irrelevant. But the 0(1, 2)/0(2) x Z2 model is different. It has no 
phase transition, and important quantities do depend sensitively on the choice of action. 

5 The lattice currents 

In the continuum formalism of paper I the form of the conserved currents was deduced from 
the group invariance of the continuum Lagrangian. The lattice analogs of the currents may 
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be deduced similarly, from the group invariance of the lattice Lagrangian 



(5.1) 



I pos 



Here the summation is over the "positive" values of A, namely (1, 0, 0, 0), (0, 1, 0, 0), (0, 0, 1, 0) 
and (0, 0, 0, 1), and the factor rj^ assumes the values —1, 1, 1, 1 correspondingly. Introduction 
of f]^ and the overall minus sign means that we are working for the moment on a Minkowski 
lattice. The reason for this is to allow the lattice field equations to have nontrivial solutions 
so that there will be something nontrivial to be conserved! 

Reference to eq. (1.4.1), with the index a replaced by the pair ab, shows that group 
invariance of Let may be stated in the form 



= 6La 
1 



i pos 



■Q\b (M + 



^a+A 



Q 



ab 



^(x+X 



Combining this with the lattice field equat 
dS 5 V- r 



;ions 







'Jo. 



-1^ E ^A 



a' 



A' v^a 



A pos L 

and remembering that the 5^"^ are arbitrary, one gets 



o- ^ vx 

^ A, 



k pos 



>a+A 



= E 

A pes 



^a+A 

3ab 



■Q'-b (0a+A) - 



Jab^ (a + ^a) - jab-^ (a - 



where 



Jab^ (a + ^a' 



■Q\b (0a) 



■Q'ab {K+x) ■ 



(5.2) 



(5.3) 



(5.4) 



(5.5) 

'^a+A 

Equation ( p.4| ) is the lattice conservation law and expression ( |5.5| ) is the lattice current 
vector. 
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Passage to the Euclidean lattice is effected by setting 77^ = 1 for all A and writing 



JabX ['^+2 ) ^ ['t'c+X) ■ (5-6) 

In paper I an alternative and simpler form (1.4.11 or 1.4.17) is given for the currents, based 
on the action (1.2.5). A Lagrange multiplier can also be introduced here, with the result 
that expression (5.6) gets replaced by its equivalent 



2 ^ 2 "'^g+A 

Making use of the explicit forms (1.3.7) and (1.4.15) one finds 



smn lA I </>Q,, (pct+X) 



and, in particular. 



1,\ _ n ^{^a^'Pa+X 



The lattice Fourier transforms of these currents, namely 

5„.aW = 'V-'i:^„Af« + iA)-"'*""'^'. (510) 



are to be used in the lattice versions of eq. (1.10.19): 

Jo.xi^Kx'(^'))^,Pn5, ^' ""^^^^ ""^^^ S,^-^' , A ^ A'. (5.11) 
where 

PR = f-iW. (5.12) 

Note: Only not fiR, can be obtained from the computer. If hr is to remain finite as 
a — > then must vanish in the continuum hmit. 

6 The small- /3 approximation 



When the lattice action is chosen in the form ( |4.2| ) the lattice version of eq. (1.10.5) may be 
written 

j A[(j)] exp Efy^ x 0q;+a)] 11 ^inh sa dsa dOoc 

{A[(P\) = ■ ^ (6.1) 

/ exp Ea,X (<^ck' ^ct+x)] 11 ^inh sa dsa. dOa. 
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where (cf. eq. (1.3.8)) 



cosh A 



'a^^a+X) = cosh cosh s^^;^-sinhs^sinhs^^;^ cos (^Q,-^^^;^) . (6.2) 

When /3 is very small the integrals in (|6.1|) receive important contributions from a vast range 
of fields. For most of these fields A fluctuates over very large values. If the vacuum is fixed, 
with the lattice averages of (p^ and 0^ held equal to zero, the sa too fluctuate over large 
values, and eq. (|6.2|) may be accurately replaced by the approximation 



2*^ 



i exp (s^ + s^^x) [l - cos (^a ' ^oc+X 



whence 



(6.3) 



(6.4) 



The corresponding approximation for the expression standing under the integral sign of 
the denominator in eq. (|6.1| ) is 



exp 



1 



1 



Y\ sinh Sex dscx dOcx. 




The first term in the final exponent represents the effect of the vast area of configuration 
space that must be probed by the functional integrals as /? — 0. The linearity (in the sa) 
of this term suggests a shift of zero point. Indeed, under the replacement 



expression ( |6.5| ) becomes 

Ar4 



,1/32/3 



1 UOcx 



( - ^ E I (^a + ^cx+xf + \ In sin' \ (O^ - 9^^^ 
^ cx,X ^ 



+ 2s/p {u^ + u^^x) In sin' I {9^ - 9^^^ 



(6.6) 



+ /3 



n 2 



Insin^- 



^(x+X 



l[dsad9a. (6.7) 



cx 
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In the ratio ( |6.1j ) the factor in front of the above expression is irrelevant, and if the terms in 
a/5 and (3 in the exponent are discarded, the action takes the effective hmiting form 



4 , 

OL,X 



(6.8) 



which is independent of (3\ Even more remarkable, the variables ua are governed by Gaussian 
statistics and decouple from the variables 9ct- 

The above simple steps appear to give us direct access to the continuum limit. We must, 
of course, check this, through use of the exact lattice action. The computer results show, 
in fact, that the approximations embodied in eqs. (|6.4| ) and (|6.8|) give a good description of 
the 0(l,2)/0(2) X Z2 model in the small-/? limit. For the rest of this section and the next 
section, therefore, we investigate the implications of the lattice action (|6.8|) . 

The cu-part of S^j^e is easily expressed in terms of Fourier components. Proceeding exactly 
as in the derivation of eq. ( p.lOp one finds 

lE K+^c.+A)' = ^E^'(k)l-kr (6.9) 

where 



- . I irk Til vrm Txn\ , . 

K{k) = hcos— ,2cos-,2cos— ,2cos— j , (6.10) 

k\k) = 16-K\k). (6.11) 
From this it follows that 

(cDkcDkO = 5k,-k'/^'(k). (6.12) 

Just as in the case of the massless free field there is a zero mode, but here it corresponds 
to k = ks where k^ = Qa^, |A^, |A^, |A^) rather than to k = 0. On the lattice this mode 
has the form 

where 

|a| = |a| + + I7I + a = (a,/5,7,5). (6.14) 

Because of its "spikiness" expression ( |6.13| ) will be called the spike function. Any multiple 



of the spike function can be added to ua without changing the value of the action S^j^e- 

The spike mode must be suppressed (i.e., the sum Z]q;(~1)''^''-^q: must always be assumed 
to vanish) for two reasons: 
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1. The spike function does not respect toroidal periodicity when is odd and hence does 
not exist as a zero mode in that case. It is a spurious lattice artifact. 

2- If Z]q:(~1)''^'^^q; were allowed to drift arbitrarily far from zero then half of the Sa of 
eq. ( |6.6|) could become small enough to render invalid the approximation ( |6.4| ) which 
led to the introduction of the ujcy. in the first place. 

Suppression of the spike mode does not imply that spikiness itself is removed from the 
theory. Indeed it shows up in the 2-point function {ucxUJa.')- Using double primes on 
summation signs to indicate that terms with k = kg are to be omitted, and introducing the 
abbreviation 

k = k- k,(modAr), k{k) = K(k) , (6.15) 

one finds 

{uJcx ujcx') = Y." {^^ ^k') e-(^-^/^)(^-"+'^'-"') = n{n - a') (6.16) 



k,k' 



where 



-(27ri/7V)k-a 

k ^^(k) 



g -(27r»/JV)(k+k 



N-'Y .-../.X = (-l)l"l$(a) , (6.17) 



k KH\. 



$(q:) being the free-field 2-point function defined in the Appendix. Since $(q;) is a smooth 
function VL{cx) is a spiky function. Its spikiness indicates that there are correlations between 
the cj's at neighboring sites analogous to spin correlations in antiferromagnetic matter. 



Equation ( |6.17| ) is strictly valid only when N is even. When N is odd it should be 
replaced by 

fi(a) = (-l)l"l$i/2(a) (A^odd) (6.18) 

where $1/2 is defined by the sum ( [A.2|) in the Appendix but with the components of k 
running over half-odd-integral values. VL{ol) is seen to be spiky for both even and odd A^. 
On the computer it is convenient always to work with A^ even. One can then separate the 
lattice into even and odd sites, and this permits efficient computer programming through 
"vectorization." 
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Equations ( |6.6|) and ( |6.16|) lead to an easy derivation of (<^°) in the small-/3 limit, pat- 



terned on the corresponding derivation (|3.14| ) for the 0(1, l)/2'2 x model: 



The last line follows from 



cosh ^ ^cosh + sinh ^ ^sinh -^^ 

1 1/16/3 ^ /^2\" ^ ip[|+f^{0)]/2/3 

2' (2/5)^^^! / 2' 

(6.19) 



^](0) = $(0) = F(Ar) — ^.1549... (6.20) 



see eqs. (|6.17| ) and (|A.5|) ). (0°) here, like (0°) for the 0(1, 1)/Z2 x Z2 model, has an essential 



singularity as /5 — ^ 0, which makes it difficult to evaluate on the computer. It is, in fact, 
doubly hard to compute for the following reason: In the small-/3 limit a histogram of the 



accumulated values that s assumes over the lattice has a Gaussian form of width ^^VL{^)j(i 
peaked at 1/16/3 = .0625//5. (See section 8, Figs. 2 and 3, for numerical confirmation of 
this.) But the most important contributions to (0°) come from values of s near .14//5, which 
is way out in the tail of the Gaussian distribution. The computed value of (0°) is extremely 
sensitive to the tail of the histogram, and reliable values for it simply cannot be obtained. 
Consequently, eqs. (1.8.20) to (1.8.23), which are based on the fields 0", do not provide a 
useful route to the renormalized Planck mass or, on the computer, to the constant (5^ of 
eq. (^ ). 

7 The fields ui and the ^-part of S^^q 

We have found that the Riemann normal variables a* of eq. (1.6.15) do yield reliable 2-point 
functions on the computer. In the small-/? limit the renormalization constant Z„ of 
eq. (1.10.15) becomes 

I I ( ^ s coshs \ \ 1 ^ , 



and the renormalized fields (1.10.17) are given by 
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The renormalized 2-point functions therefore have the small-/? limits 
[(^ROc^^Ra') = (32/3)^ {sa Scx') (cos 6*0; cos 6*0;/) 



(7.2) 



(32/?) 



1 UoL 

+ 



1 



16/5 ' y/P J \16P ' y/P 

4 (cos 6*0; cos 6'q;') 



(cos^^q: cos^^q;') 



and, similarly, 



<^Ra<^Ra'} ^ 4 (cos^a sin^^') 



2 2 



Ra^fio:' / ^ 4 (sin ^q, sin 6 a') . 



(7.3) 

(7.4) 
(7.5) 



Owing to the 0(2) invariance of the vacuum one expects expression (7.4) to vanish and 
expressions ( [7.3| ) and ( |7.5| ) to be equal, so that 

\ .^..f^f^. (7.6) 



where 



9(q: — a') = 4 (cos 6'q; cos 6'q;') = 4 (sin 6'q; sin 6'q;') 

= 2 ((cos ^Q cos 6'q' + sin 6'q sin 6'q')) = 2 (cos (6'q; — 6'q')) . (7.7) 

The question of the existence of a continuum limit for the 0(l,2)/0(2) x Z2 model is seen 
to be determined entirely by the 6'-part of the action Suj^e- The reader should note that the 
above limiting forms for the renormalized 2-point functions are all independent of /3. 
The 6'-part of S^^e gives rise to the following unusual average: 



de 



-1/4 



/^^Oq.a 

n/ dOoc. 



pos 



sm \ (e^ - e^^)^ 



-1/4 



(7i 



(7.9) 



Note that the integrands become singular when the cross-link differences = 6^ — 6fy^^\ 
vanish. One therefore expects histograms of the A^'s to be strongly peaked around zero. 
This is confirmed by the computer results. (See Fig. 1.) The singularity of the integrands 
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stems from the fact that the action S'^g is unbounded from below (see eq. (|6.8|)). The 
question immediately arises why the Monte Carlo procedure works when S"^ e is used. It is 
shown in the Appendix that the singularities to which S^^^o gives rise are just soft enough for 
the functional integrals to exist. 

Another unusual feature of the average ( |7.8| ) is that its weight, or "measure," in addition 
to being invariant under the obvious 0(2) transformations 

0(x^ect + constant, Oa ^ -Oct , (7.10) 

is invariant as well under a set of nonlinear transformations having the infinitesimal form 

69a. = ^1 cos 9a + €2 sin 9(x . (7-11) 

The existence of the latter invariance, which is demonstrated in the Appendix, derives di- 
rectly from the 0(1,2) invariance of the measure of the average ( p.l|) . The transformation 
( [r.ll| ) is easily recognized as an expression of the well known "headlight effect" associated 
with Lorentz boosts. 

When using the "exact" action ( |4.2| ) one must program the computer to execute 0(1,2) 
boosts at the end of each Monte Carlo sweep in order to keep the lattice averages of 0^ 
and 0^ equal to zero and hence to keep the vacuum state from drifting. When using the 
6'-part of the action S^^e one must, for the same reason, program the computer to execute 
transformations of the form ( [7.11| ) after each sweep. In this case the vacuum is held fixed if 
the constraints 

COS = 0, sin 9a = , (7.12) 

a. a. 

which express the 0(2) invariance of the vacuum, are maintained. At the end of each sweep 
the 9a will be fairly randomly distributed from to 27r (if the lattice is big enough) and will 
approximately satisfy the equations 

A^-^^cos^^a = A^-^5^sin2^„ = 1, A^-^ J] sin cos = . (7.13) 
Ct a ^ cx 

Moreover, eqs. (7.12) will be violated only slightly. Defining 

AT-^^cos^a = 0<1, A^-^5^sin^a = 5 < 1, (7.14) 

Ot Ct 

one easily sees that the constraints will, to good accuracy, be reestablished by the transfor- 
mation (|7.11|) if ei and €2 are chosen to be 

ei = -2S, 62 = 20. (7.15) 
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This choice has been found to work well on the computer. 

We close this section by noting that, in the small-/? limit, not only the a 2-point functions 
( [7.6| ) but also the current 2-point functions ( |5.11|) and (5.12) are determined entirely by the 
6'-part of the action S^^g. Referring to eqs. (|5^ ), ( |6.4| ) and ( |6.6| ), one easily obtains the 
limiting forms 



■ f^^lx^ 1 ^os^^+A^cos^ 1 sin|(VA + ^ay 



. 1 sing sin 1 ^o^\[^oc+X + ^ocj . 



2 7 /3^o 16 sin^ \ {e^^^ - 6^) 8 sin \ {O^^^ - 6^ 

which are independent of the ujcx,. 

8 Computer results 

If one is to use the limiting action 5*^^ e with confidence, one must check how well the results 



obtained with the full action ( [4.2| ) tend, when (3 is small, toward those obtained with S'^j^g. 
This immediately raises the question: How small is smalll 

Consider eq. ([7.31). If [3 is not set equal to zero this equation implies 



<^RCt(^R(x') = [4 + (32)2/5r](a - ct')\ {cos 9 a cos 9 a>) , (8.1) 

and one sees that the second term in the brackets cannot be ignored unless 

/? < 4(32)-7r](0) = .025 . (8.2) 

Such small values for (3 have the potential of severely taxing the ability of the computer. In 
early runs attempts were made to compute (^(f'Ra'pRcy.'} at /3 ~ .002, without success. The 
cause of the failure was the same as that which makes the computation of (0°) impractical, 
as explained at the end of section 6. On the computer the difficulty can masquerade as 
critical slowing down, causing large autocorrelation intervals, drifting averages, etc. But in 
fact critical slowing down has nothing to do with it. Computations of (cT/jaCTijcK') have 
turned out to generate autocorrelation intervals of less than 20 sweeps, even for /3's as small 
as 10"'^. Moreover, the results have been found to be both reliable and accurate. 
A simple measure of the validity of the small-/? approximation is the quantity 

/|A'-A|\ 
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where A is the exact cross-hnk geodetic distance, A' is the approximation to it given by 
eq. (p^), and the average is over hnks and sweeps. At (3 = 10~^ expression ( |8.3| ) is found to 
be aheady less than 10~^, so the small-/? approximation is indeed a good one. 

More detailed measures of the quality of the small-/? approximation are provided by 
histograms of the scx and of the cross-link differences A6 = d(y. — (^cx.+X- Figure 1 shows 
the A^ histograms for both the exact lattice action (^4.2|) , at /? = 10"'^ and = 10, and 
the ^-part of the approximate action Su)^-, at = 10. The two curves, which were compiled 
by accumulating all cross-link differences for 60,000 sweeps, are indistinguishable. Figures 2 
and 3 show the sql histograms for the exact action, at /? = 10~^ and 10~^ with = 10, 



superimposed upon the theoretical Gaussian curve derived from eq. (|6.6| ) and the cj-part of 
the action Si^^q. Again the agreement is excellent. 

It should be stated that during the runs at /? < 10~^ the spike mode (see section 6) 
began to appear spontaneously and randomly, in the fields generated by the Monte Carlo 
procedure. At these values of (3 this spurious mode was forcibly removed from the field at 
the end of each sweep after the vacuum-fixing boost was applied. As remarked in section 6, 
suppression of the spike mode does not remove spikiness itself from the theory. (See the 
discussion of the a 2-point function below.) 

It should also be stated that during the small-/? runs the vacuum-fixing boosts applied 
at the end of each sweep were chosen so as to enforce the conditions ( [7. 121) rather than the 
conditions 

E<^a = 0' E</'a=0- (8-4) 



a a. 

273±33 



Because of the huge values assumed by the (~ lo^'i'iio -vyj^en j3 = 10^^ and ~ 10 
when jS = 10~^) enforcement of ( ^.4|) is impractical at these /?'s. We could have imposed the 
alternative conditions 

E^« = 0' E^a = 0, (8.5) 
a a 

but we felt that, to make comparisons, it was best to uniformly impose conditions ( |7.12| ), 
which are the only ones available when the action S^^^e is used. It turns out that enforcement 
of (|7.12|) causes conditions (|8.5|) also to be satisfied to better than .03% accuracy provided 
A^ is big enough (> 10) so that a "snapshot" of the lattice after each boost provides a fair 
sampling of field values. 

The results displayed in Figs. 1, 2 and 3 give us confidence in proceeding to the com- 
putation of Pr. We discuss first the determination of based on the lattice version of 
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eq. (1.10.16): 




Rk"^ Rk' 



(8.6) 



Here a^j^ is the lattice Fourier transform of the field c^^cx — ^a^'^a.- Note that the determi- 
nation of i3r requires a computation of the renormalization constant Z^r (eqs. (1.10.15) and 
( [TTD ) as well as the a 2-point function. According to eq. (|8.6|) Pr is given by the intercept 
(at k = 0) of the function 2/K'^(k) (5")jk^i?k*)- Figure 4 shows plots of this function for 
= 10 and various values of (3, as well as of the corresponding function for the action S*^ g 
(eqs. ( [7.6|) and (|8.7| )). At the smaller values of (3 and at low momenta the curves fall right on 
top of one another. The differences between the curves at higher momenta is explained by 
the spikiness of the function (c"1jq;C"/jq:') arising from Q{a — ck') in eq. ( ^.1| ). The spikiness 
is a short-wavelength phenomenon that is picked up by the Fourier transforms. It goes away 
as /? ^ 0. 

Figure 5 shows a plot of (3^ vs. (3 as determined from the intercepts of the curves shown 
in Fig. 4. It is on this plot that we base our conclusion that the quantized 0(1, 2)/0(2) x Z2 
sigma model has no continuum limit in four dimensions. The key fact is that (3^ vanishes 
nowhere, not even at f3 = 0. In view of the definition of (3^ (eq. (5.12)) this implies that the 
renormalized Planck mass /ir must become infinite in the continuum limit a — > 0. 

Figure 6 shows a plot of vs. A^ for the action S*^ i.e., for /5 = 0. It is evident that Pr 
tends to a constant value as A^ increases past 10. Its value on an infinite lattice is therefore 
effectively known, and one can pass to the continuum limit a without worrying that 
the size of the "universe" (the lattice 4-torus) is simultaneously collapsing to zero. The 
nonvanishing of Pr, and the blowup of fiR, is not an artifact of a vanishing universe. 

If a similar result (infinite (ir) is eventually found in conventional quantum gravity we 
shall be forced to conclude that the theory is not viable, because the classical theory of 
Einstein cannot emerge from it as a low-energy limit. In the case of the 0(l,2)/0(2) x Z2 
model the blowup of fiR implies that there can be no low-energy "classical" regime in which 
the physics is described by an action functional having the form of the first term on the 
right side of eq. (1.8.18). It is in this sense that we say that the quantized 0(1, 2)/0(2) x Z2 
model has no continuum limit. 

It may occasionally be useful to adopt an alternative viewpoint by saying that the contin- 
uum limit exists but is trivial. However, this viewpoint can be misleading on two accounts: 

1. It reinforces the current folklore which claims that all scalar field theories in four 
dimensions are trivial. The claim is suspect because it is just folklore. No valid 
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argument exists that would have given us knowledge a priori that /i/j must blow up 
in the continuum limit. Each scalar theory must be investigated on its own. For 
example, despite the amount of research that has been done on the 0(3)/0(2) model, 
it is still unknown whether a continuum limit exists for it when the action ( |4.2| ) is 
used. If does vanish anywhere for this model one presumes that it will do so at the 
phase-transition point /? = Pc, but no one knows. 

2. While it is true that the 0(l,2)/0(2) x Z2 model becomes trivial in the continuum 
limit, in the sense that the renormalized "coupling constant" /i^^ vanishes, it is also 
true that in this limit the model is no longer 0(l,2)/0(2) x Z2. The configuration 
space becomes flat and the invariance group consists of the isometries of the plane, not 
the elements of 0(1,2). Therefore, although the continuum limit "exists," the model 
does not. 

We remark that the continuum limit may properly be regarded as located at (3 = 0, not 
only because the naive continuum limit is located there (with the bare Planck mass /x held 
at a fixed finite value) but also because the lattice specific heat diverges there, making (3 = 
a critical point of the theory. This may be seen by introducing the "partition function" 



Z(/3) 



fi[(p] d(j) 



•7) 



and noting that eq. ( p. 201) is equivalent to 

1 cPZ 



CL = N- 



1 

Z2 



'dZ' 

dd 



In the small-/? limit the partition function is dominated by the first factor in expression (|6] 
whence 



This immediately yields 



cl 



const, xe 



const. x(3 



7V4/32/3 



(8.9) 



^.10) 



which diverges at /3 = 0. The critical exponent, —3, in eq. ( |8.1CI|) has been confirmed by 
direct computation of ( |3.20| ) for the exact action ( |4.2| ). Note that this critical exponent 
differs from that of the massless free field (eq. ( |3.21| )). 

The "triviality" of the continuum limit of the 0(1, 2)/0(2) x Z2 model, which is entirely 
due to the nonvanishing of (3r, implies nothing about the shape of the /3 = curve in Fig. 4. 
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Since this curve is plotted effectively as a function of v^j iir ■, and since [iji diverges, finite 
physical values of are all located at the left-hand intercept. It is therefore quite remarkable 
that this curve is nevertheless very nearly a straight horizontal line. Such a result was by no 
means obvious a priori. Suppose one were to present a mathematician with eq. ( [7.8|) without 
telling him where it came from. How could he possibly know in advance that the average 
([7.7|), defined by ( |7.8| ), would be very nearly the 2-point function for a pair of free fields? 

There is a consistency check on this curious result. Once one knows that the average 
( [7.7| ) is very nearly proportional to the 2-point function $(q; — a') for the free field, one may 



use the proportionality constant to define [3r aX. [3 = 0: 

Q{a - a) K I3r^<^{(x - a') . (8.11) 
But from eq. ( |7.7| ) it is obvious that 0(0) = 2. Moreover 



whence it must follow that 



$(0)=F(iV) — > .1549, (8.12) 



Pr^\f{N) .0775, (8.13) 



which agrees very well with Fig. 6. 

It is useful to sketch the procedure that we should have had to adopt if things had turned 
out differently and I3r were to vanish for some value of (3. To begin with we should have 
tried to match the a 2-point function of the continuum theory to an Ansatz of the form 

or some variant thereof. Here a and A are dimensionless fitting parameters independent 
of /iij. Equation ( ^.14|) would be a universal formula, with fj,ji merely setting the scale. 



The lattice analog of (|8.14|) is 



^rA') = r . (8.15) 

This formula becomes meaningless when (3r vanishes, and one might therefore be tempted 
to conclude that, since the average (|7.7|) is independent of /3, it was obvious, as soon as 
eqs. ( [7.6|) and (|7.7|) were written down, that [3r cannot vanish. But this would be too hasty 
a conclusion. Every finite lattice has associated with it an "infrared" mass 

y,' = l/L = l/Na. (8.16) 
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The corresponding lattice quantity is 



11 a = N 



(8.17) 



which, because it remains finite as a goes to zero, could quite possibly mask jS^ in the limit. 
If the real (Sr had turned out to vanish then one would have found that the effective I3r 
defined by eq. (8.6) was proportional to A^"^. One might have regarded any hope for such a 



result as foredoomed, in view of the insensitivity of most averages to A^, once is bigger than 
6 or 8. But eq. (|7.8|) defines an unusual kind of average, and we had no advance knowledge 
of its properties. 

We turn finally to the determination of (3r from the current-current average ( ^.11| ). Be- 
cause of large scatter in the individual values of jQiX^^jQ^xi.^)* "we found that it was not 
possible to use eq. (|5.11 ) at values of (3 below 10^^. But down to that value the agreement 
between the values of (3r obtained from the a 2-point function and those obtained from 
( p.ll| ) is good. Indeed, Fig. 7 shows that the agreement is remarkably good. Because the 
determination of (3r from ( p.ll| ) avoids computation of the renormalization constant Z^j, this 
good agreement is far from trivial. It is a convincing demonstration of the consistency of 
the formal arguments presented in paper I. It makes one believe that there is something to 
quantum field theory after all! 

In closing we note yet another consistency check on the formalism. Equations (1.9.2), 
(1.9.3), (1.9.4) and (1.9.7) together imply 



(vac| T (a\x)jQj^{x)) |vac) 



P 



Replacing a* by cr}j, passing to Euclidean space and taking Fourier transforms, we get 



(^k(p)Jo/(p) 
The lattice version of this equation yields 

(4kJoa(k)* 



i5i~5{p + p') 



X ■ K(k) 



(8.19) 



^.20) 



which is a universal function, the same for all values of /3. Equation (|8.2CI|) has been checked 
at = 4 for all values of k and A and for many values of (3. We have found it to be satisfied 
to 5% accuracy for values of /3 down to .05, and to 20% accuracy at (3 = .02 . At lower values 



of P large scatter makes checking impossible. Although eq. ( p.20|) gives no information about 
the continuum limit, it is nevertheless a remarkable equation. Its left-hand side depends on 
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Zf, and /3; its right-hand side depends on neither. Its "experimental" verification gives us 
one more reason to beheve that we have done the right things. 
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Appendix 

Estimating F{N) 

To estimate F{N) use the approximation ( p.lSI) . The exact equation for the 2-point function 



IS 



where (cf. eq. ( |2.14| )) 



Evidently 



and, from eq. (2.15), 



, ^-{2TTi/N)k-CX 

$(0) = F{N), 



1 



.025 ... for all A . 



Now use the fact that $ satisfies a simple equation involving the lattice Laplacian: 
1 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



[-"^{cy. + A) + 2$(ck) - $(a - A)] 



A 

1 



A k L 



(A.5) 



in which eq. (|2.13|) is used in passing to the last line. Setting a = in ( |A.5D and remembering 
that A ranges over eight different values, get 
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Invariance of the measure of the average ( |V.»|) 

The measure may be written in the form 
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Setting 9'(y^ = 6(x+S6a. one finds for the change in d6cx under the infinitesimal transformation 



6d9c = dd'^ - ddct = [{d9'c,/d9a) - 1] dOa = {d69a/d9a) dOcc 
= (-ei sin 6a. + £2 cos 6(x) dd^ 



and hence 
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the trigonometric identities 

cos X — COS y = — 2 sin I (a; — sin ^ (a; + y) 
sin X — sin y = 2 sin i (x — y) cos | (x + 

being used in passing to the final form. The complementary identities 

cos X + COS y = 2 cos ^{x — y) cos ^{x + y) 
sin X + sin ?/ = 2 cos | (x — sin | (x + 
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permit expression ( |A.1CI| ) to be further simplified, and one sees tliat, under ( [7.11| ), tlie measure 
( |A.7] ) gets multiplied by 

^ < -ei sin 9a + £2 cos Oa - j^Yl (^"^ + ^cx+x) + ^2 (cos 9^^^ + cos 6'q,_^;^ 
« I ^ A 

(A.13) 

But this vanishes by virtue of the fact that A ranges over 8 values per site. It follows that 
the measure is invariant under ( 7.11| ). 



Proof of the existence of the functional integrals in (\7.i^ ) 

The integrands become singular when two or more 6*0; 's approach simultaneously a common 
value 9. Consider therefore a rectangular subvolume of the lattice containing JKLM sites, 
where J, K, L, M are integers between 1 and A^, at least one of which is greater than 1. 
The number of links in this subvolume is AJKLM - JKL - J KM - JLM - KLM. Replace 
the integration variables 9oc attached to the JKLM sites by ^, A9 and JKLM — 2 angular 
variables in {JKLM — l)-dimensional Euclidean space. The variable A9 is the maximum 
distance of any of the 6*0; 's from the common value 9 and is the critical variable for de- 
termining the convergence of the associated subintegration. At small values of A6' the A^^ 
subintegration goes like 

= / (A^)K-^^^+^^^+-^^A^+™)-2 d(A9) . (A. 14) 

Jo 

But JKL + J KM + JLM + KLM > 7. Therefore this subintegration is no more singular 
than Jq{A9)~^^'^ d{A9), which converges. Since all the subintegrations converge, the integral 
for the whole lattice converges. 

The convergence is no accident. All the ^a's on the whole lattice can be brought arbi- 
trarily close to one another by imposing a sufficiently powerful boost (the headlight effect). 
But the integration measure has been shown above to be boost invariant. 

It is a curious fact that the measure ([A.7|) is the same as the measure that appears in 



the computation of the string amplitude for A^^ open-string tachyons on the boundary of 
a disk. It is a consequence of the 0(1,2) symmetry of both systems. The above proof of 
convergence is the same as that used in string theory. 
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Figure Captions 

Figure 1. Histograms of the cross-link differences A^^ = 9^ — for the exact lattice 

action ( ^?2|) , at /3 = 10""^ and = 10, and for the 6'-part of the approximate action S^^g 
(eq. (|6.8|) ) at = 10. All cross-link differences for 60,000 sweeps have been compiled into 
101 bins. The normalized distributions for the two cases are indistinguishable. 

Figure 2. Histogram of the sa at /3 = 10^^ and A^ = 10, superimposed upon the theoretical 
Gaussian distribution derived from eq. ( |6.6|) and the cj-part of the approximate action 8^0,0 ■ 

Figure 3. Histogram of the sct at /5 = 10~^ and A^ = 10, superimposed upon the theoretical 
Gaussian distribution derived from eq. ( |6.6|) and the cj-part of the approximate action S^j^e- 

Figure 4. Plots of i^Rk^Kk*) vs. K^Pr for various values of P, a.t N = 10. Pr is 

determined for each curve from its intercept at = as obtained from a linear fit to the 
low-momentum values. The curves are plotted as functions of K'^/Pr because since K'^/Pr 
becomes /i/j^ in the continuum limit, this is the only scale- invariant way to compare them. 
Although the curves are not perfectly smooth they show remarkably little dependence on 
the detailed values of the components of k. Their essential dependence on alone shows 
that 0(4) invariance is maintained in the continuum limit. (See the discussion following 
eq. (|2.15|) .) The error bars (statistical plus autocorrelation) range from .1% at the high- 
momentum end to 1% at the low-momentum end. Only the latter are shown on the figure. 
For P < 10^^ the curves are seen to fall on top of one another at low momenta. The difference 
between the curves at higher momenta is explained by the spikiness of (d^Rk^Rk*) arising from 



il{a — a') in eq. (|8.1|) . The spikiness is a short- wavelength phenomenon that is picked up by 



the Fourier transform. It goes away as /? — >• 0. The /3 = curve itself is obtained by using 
the 6'-part of the action S^^^q. 

Figure 5. A plot of Pr vs. P as determined from the intercepts of the curves shown in 
Fig. 4. The error bars, which are less than a tenth of a percent, are almost too small to be 
seen. At higher values of P the model tends toward a free-field model, and Pr p. It is 
clear that there is no value of P for which Pr vanishes. 

Figure 6. A plot of Pr vs. A^ when P = 0, obtained using the 6'-part of the action S^j^- The 
asymptotic value of Pr, as N 00, is approximately .078 . 

Figure 7. Plots of Pr vs. P at N = 10. The black circles represent the values obtained from 
the a 2-point function. The triangles represent the values obtained from current-current 
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averages at low momenta. Error bars are shown only for the triangles. Error bars for the 
circles are too small to be seen. 
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